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Abstract 

The reaction of cyclobutane-l,2-dione with hydroxide was studied by a variety of ab initio (MP2, SCS-MP2, CCSD(T), CEP A/1) 
and density functional (M06-2X) methods. Three possible reaction paths of the initially formed tetrahedral adduct leading to either 
1-hydroxycyclopropane-l-carboxylate (benzilic acid type rearrangement, path A), a-oxobutanoate (path B) or y-oxobutanoate (path 
C) were considered. Although the latter two products show similar or even more negative Gibbs free energies of reaction than 
calculated for the benzilic acid type rearrangement, the Gibbs free energies of activation are substantially higher. According to the 
calculations, the only feasible reaction appears to be the formation of 1 -hydroxy cyclopropane- 1-carboxylate, which is corroborated 
by previous experimental observations. 

Introduction 

Addition of nucleophiles, e.g., OH~, to 1,2-dicarbonyl com- 
pounds leads to the formation of relatively stable tetrahedral 
adducts (Scheme 1) [1]. These adducts further react either by (i) 
fission of the R 2 -C bond and migration of R 2 (benzil-benzilic 
acid rearrangement, path A); (ii) fission of the R 2 -C bond 
without migration of R 2 resulting in a-oxocarboxylic acids 
(path B); or (iii) fission of the carbonyl-C-sp 3 -C with forma- 
tion of an aldehyde and carboxylic acid (path C). In the case of 
benzils, depending on the substituents on the aryl rings, all three 
types of reactions have been observed [2]. 
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Scheme 1: Reactions of 1,2-dicarbonyl compounds with base. 
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The benzilic acid rearrangement of cyclic 1,2-diones [3,4] leads 
to ring contraction, e.g., the rearrangement of cyclobutane-1,2- 
dione (1) to 1-hydroxycyclopropanecarboxylate (2) [5-7]. In 
contrast, cyclobut-3-ene- 1,2-diones 3 react to 2-oxobut-3- 
enoates 4 (at least formally according to path B) [8], whereas 
benzocyclobutene- 1,2-diones 5 lead to 2-formylbenzoates 6 
(path C, Scheme 2) [9]. 




Scheme 2: Reactions of cyclic 1,2-diones with base 



In view of the importance of the benzil-benzilic acid rearrange- 
ment in organic chemistry, several computational studies 
concerning this rearrangement [10,11] or related reactions 
(Favorskii rearrangement [12-14], halolactonisation [15]) have 
been published. To the best of our knowledge no attempt has 
been made so far to consider the additional pathways B and C in 
these reactions. Here we present a detailed computational 
study (DFT and ab initio) of the base-catalyzed reactions of 
cyclobutane-l,2-dione (1) taking into account all three possible 
pathways. 

Results and Discussion 

The various transition states, intermediates and products 
initially considered for the three reaction paths A, B, and C are 
depicted in Scheme 3. It turned out that not all of the structures 
shown in Scheme 3 could actually be located as stationary 
points on the potential-energy surface. On the other hand, some 
other stable as well as highly reactive intermediates and/or tran- 
sition states were obtained (see below). Generally, in nucleo- 
philic addition reactions to carbonyl compounds in aqueous 
solution, water not only acts as a solvent but frequently actively 
participates in the reaction, such as in water-assisted hydrolysis 
[16-20] or also in the benzil-benzilic acid rearrangement [11]. 
Therefore, cyclobutane-l,2-dione (1) hydrated by two water 
molecules, and hydroxide ion hydrated by four, i.e., 1(^0)2 
and [OH(H20)4]~ , were used as reactants. Hence for all tran- 
sition states, intermediates and products shown in Scheme 3, 
hydration by six water molecules is implied. 
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Scheme 3: Possible intermediates, transition structures, and products 
considered for the reaction of cyclobutane-1,2-dione with hydroxide 



Relative Gibbs free energies with respect to the separated reac- 
tants 1-(H20)2 + [OH(H20)4]~ including bulk aqueous solva- 
tion energies (SMD solvation model [21]) obtained by various 
computational procedures are collected in Table 1. More 
detailed results are provided in Supporting Information File 1 . 
Before discussing the individual reaction paths in detail, a com- 
parison of the computational procedures used (M06-2X [22], 
MP2 [23] and SCS-MP2 [24], the composite energy scheme 
E c = £(MP2/6-311+G(2df,2p) + £{[CCSD(T) - MP2]/6- 
31+G(d)}; and LPNO-CEPA/1 [25-27]) is made. 

Comparison of computational procedures 

Taking the CEP A/1 results as a reference, inspection of Table 1 
reveals that all other computational procedures invariably lead 
to a greater stabilization of all stationary points considered with 
respect to the separated reactants. The largest deviation from the 
CEP A/1 energies is obtained with M06-2X, the smallest with 
SCS-MP2. The corresponding mean absolute deviations and 
RMS errors are 5.9, 6.3 (M06-2X); 3.6, 3.9 (MP2); 2.5, 2.9 
(SCS-MP2); and 3.9, 4.0 (E c ). However, the predicted trends 
agree very well with that obtained with CEPA/1 [Supporting 
Information File 1 ; the corresponding correlation coefficients 
R 2 are 0.992 (M06-2X); 0.997 (MP2 and SCS-MP2); and 0.999 
(Ec)]. Inclusion of an empirical dispersion correction (DFTD3 
[28]) in the M06-2X results further lowers the relative energies 
by ca. 0.8 kcal mol -1 . 

Figure 1 summarizes the energetic as well as structural aspects 
of all three possible reaction paths in the reaction of cyclo- 
butane-l,2-dione with hydroxide anions. 
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Table 1 : Relative Gibbs free energies for all stationary points consid- 
ered 3 . 

M06-2X MP2 SCS-MP2 E c CEPA/1 



Addition of neutral nucleophiles to carbonyl compounds gener- 
ally proceeds via formation of a prereaction complex; if Gibbs 
free energies instead of pure electronic energies or enthalpies 



are considered, such complexes usually become unstable [29]. 
In the present case without inclusion of bulk solvation effects 
(Supporting Information File 1) this prereaction complex 
(ion-dipole complex IDC) is quite stable with AG = -14 to 
-19 kcal mol -1 compared to the separated reactants. However, 
especially the large calculated [SMD-M06-2X/6-31G(d)] solva- 
tion energy of [OH(H20)4]~ compared with that of the 
ion-dipole complex (AG so i v = -66 versus -54 kcal mol -1 ) 
results in an endergonic formation of this complex (AG = 
2.5-7.3 kcal mol -1 , Table 1). After passing transition state TS1 
(v = 144 i cm -1 ), formation of the tetrahedral adduct by attack 
of [OH(H20)4] - at CI leads to the first truly stable intermediate 
Intl. This intermediate is common to all further possible path- 
ways. The nucleophile involved in the formation of Intl actu- 
ally is a water molecule, since formation of the CI— 03 bond is 
accompanied by proton transfer H3 to 04 of the original 
hydroxide anion (Figure 2). An analogous concerted 
addition-proton transfer has also been calculated previously for 
the benzilic acid type rearrangements of biacetyl and benzil 
[11]. Car-Parrinello molecular dynamics simulations of the 
hydrolysis of formamide in basic solution indicated that the 
traditional view of attack by hydroxide anion rather than 
a first-solvation-shell water molecule is more likely; 
however, the more powerful electrophile methyl formate 
should react according to the general-base mechanism [30]. 
The second carbonyl group in 1 ,2-diones efficiently enhances 
the carbonyl reactivity to make attack by water the preferred 
mode of reaction. 
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a AG (1 mol L~ 1 standard state, kcal mol -1 ) with respect to separated 
reactants 1-(H 2 0) 2 + [OH(H 2 0) 4 ]~. All geometries optimized with M06- 
2X/6-31+G(d,p). For M06-2X, MP2, and SCS-MP2 single-point calcula- 
tions the 6-31 1 +G(2df,2p) basis set was used; LPNO-CEPA calcula- 
tions were done with the def2-QZVPP basis set; E c = E[MP2/6- 
31 1 +G(2df,2p)] + E{[CCSD(T) - MP2]/6-31 +G(d)}; AG S0 | V calculated 
by SMD M06-2X/6-31G(d). 
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Figure 1: CEPA-1/def2-QZVPP calculated reaction paths for the reaction of 1(H 2 0) 2 + [OH(H 2 0) 4 ] _ . 
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TS1 



Int1 




r(C1-01) 1.234 r(C1-C2) 1.560 

r(C1-03) 1.867 r(C2-02) 1.205 

r(03-H3) 1.464 r(C1-C4) 1.555 

TS2 




r(C1-01) 1.254 r(C1-C2) 1.466 
r(C1-03) 1.391 r(C2-02) 1.258 
r(C2-C4) 1.863 r(C1-C4) 2.112 




r(C1-01) 1.322 r(C1-C2) 1.564 
r(C1-03) 1.445 r(C2-02) 1.206 
r(03-H4) 0.986 r(C1-C4) 1.573 



4 




r(C1-01) 1.246 r(C1-C2) 1.514 
r(C1-03) 1.404 r(C2-02) 1.272 
r(C2-C4) 1.503 r(C1-C4) 2.585 



Figure 2: M06-2X/6-31+G(d,p) calculated structures of stationary points along the benzilic acid type rearrangement. Distances are given in 
angstroms (A). 



Benzil-benzilic acid rearrangement (Path A) 

With respect to intermediate Intl, the Gibbs free energy of acti- 
vation of the benzilic acid rearrangement (fission of the C1-C4 
and concerted formation of the C2-C4 bond) via TS2 
(v = 246/ cm -1 ) is ca. 15 kcal mol -1 . The initial product Pla 
(1-carboxycyclopropanolate) is expected to easily convert to the 
final product 1 -hydroxy cyclopropanecarboxylate (2) by a 
simple acid-base equilibrium (protonation of the alcoholate, 
deprotonation of the carboxylic acid, ApK a ~ 10). Overall, path 
A is not only strongly exergonic [AG react = -24.4 (CEP A/1) to 
-33 kcal mol _1 (M06-2X)] but has also a quite low barrier 



[AG^ = 6.9 (MP2) to 14.3 kcal mol" 1 (CEP A/ 1)], Table 1. Thus, 
this reaction should be quite feasible. Pertinent geometrical data 
are collected in Table S3 of Supporting Information File 1; 
M06-2X/6-31+G(d,p) optimized structures of relevant station- 
ary points along path A are depicted in Figure 2. 

The breaking C1-C4 bond is stretched from 1.573 A in Intl to 
2.1 12 A in TS2 (2.585 A in 2), while the newly formed bond 
C2-C4 is shortened from 2.125 A in Intl to 1.863 A in TS2 
(1.503 A in 2). The Wiberg bond indices resulting from an 
NBO analysis [31] for the C1-C4 bond are 0.92 (Intl) and 0.35 
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(TS2); those for the C2-C4 bond are 0.04 (Intl) and 0.40 
(TS2), indicating nearly equal bond breaking and formation. 
The feasibility of the carbanion [l,2]-shift in the benzilic acid 
rearrangement has been attributed to the special shape [11] of 
the LUMO of 1,2-dicarbonyl compounds. 

Path B 

Product P2 should be even more stable than 2. However, 
despite several attempts neither TS3 nor Int2 could be obtained. 
Instead, invariably TS2 or Pla were obtained. It is tempting to 
assume that the relatively close contact between C2 and C4 
preferentially results in C2-C4 bond formation rather than addi- 
tion of a proton to Int2. To address this problem, optimization 
of an extended conformation Int2' was attempted. However, 
such a structure collapsed upon geometry optimization in a 
concerted proton transfer-nucleophilic addition reaction to 
intermediate Int2a. By a simple acid-base equilibrium (alcohol- 
ate-carboxylic acid — > alcohol-carboxylate), this intermediate 
is expected to convert immediately to product P2a, i.e., the 
hydrate of product P2 (Scheme 4). 



_z-^/COOH 


/^COOH 


^xoo- 


O 


HO 0" 


HO^OH 


Int2' 


Int2a 


P2a 


Scheme 4: Reaction sequence calculated for an extended con- 
formation of Int2. 



While all attempts to locate transition state TS3 as well as those 
ones leading to either Int2' or Int2a were unsuccessful, a path 
[TS3a (v = 410/ cm -1 , Scheme 3 and Figure 3] directly 
connecting Pla instead of Intl with P2, could be obtained. 
Hence, path B actually does not start off from Intl but diverges 
at the initially formed product Pla of path A. 

Compared with TS2, transition state TS3a (Figure 3) is charac- 
terized by significantly longer C1-C4 and C2-C4 bonds, 2.975 
and 2.259 A in TS3a versus 2.112 and 1.863 A in TS2. Simi- 
larly, the C2-C3-C4 angle and the C1-C2-C3-C4 dihedral 
angle are much larger in TS3a (96.7° and 76.0°) than those 
found in TS2 (76.6° and 54.8°), Table S3 in Supporting Infor- 
mation File 1 . The main "movement" in TS3a is transfer of a 
proton from 04 to C4; consequently the imaginary frequency of 
TS3a is larger (v = 410/ cm" 1 ) than that of TS2 (v = 246/ cm" 1 ) 
with C-C bond formation as the associated mode. 

PathC 

Similar to path B, all attempts to locate the initially proposed 
intermediate Int3 were unsuccessful. Instead a more compli- 



TS3a 




r(04-H4) 1.274 r(C2-C4) 2.259 
r(C4-H4) 1.554 a(C2-C3-C4) 96.7 
t(C1-C2-G3'C4) 76.0 



Figure 3: Calculated structure of transition state TS3a. Distances are 
given in angstrom (A), angles in degrees. 

cated pathway involving a high-energy bicyclic intermediate 
Int4 was found (Figure 4). 

Furthermore, this intermediate did not react to product P3 but 
instead through a concerted ring opening and proton transfer to 
P3a. The anticipated product P3 of path C is a y-oxocarboxy- 
late. Such y-oxocarboxylic acids or carboxylates are prone to 
ring-chain tautomerism [32,33]. Product P3a essentially is the 
ion-dipole complex between the ring tautomer of neutral 
y-oxobutanoic acid with [OH(H20)5]~ , and hence, its forma- 
tion is completely reasonable. The actual pathway C obtained 
by the calculations is indicated in Scheme 5. 

In fact, Int4 is barely stable if Gibbs free energies and bulk 
aqueous solvation are taken into account (Table 1 and 
Supporting Information File 1) and should collapse more or less 
barrierless to the tetrahedral adduct Intl. In contrast, 
rearrangement to P3a involves a substantially higher barrier 
(40-45 kcal mol -1 with respect to Intl) than path A (ca. 
15 kcal mol" 1 with respect to Intl). TS4, Int4 and TS5 can be 
considered as bicyclic structures consisting of a 3-membered 
oxirane and a 4-membered cyclobutane ring. The two rings are 
inclined to each other as indicated by the angle a measured 
between the midpoints of the C3-C4 and C1-C2 bonds and the 
oxygen atom 01 [a = 103° (TS4), 106° (Int4), and 112° 
(TS5)]. The product P3a of path C has a largely planar five- 
membered ring structure (a = 165°). In both TS4 and Int4 the 
C1-C2 distance (1.48 A) is in the range of C-C single bonds, 
while in TS5 this bond is significantly stretched (r = 1.853 A; 
in P3a this distance is r = 2.3 17 A). In contrast the 01 -C2 bond 
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Int4 




r(C1-01) 1.429 r(C1-C2) 1.483 
r(C2-01) 1.504 r(C2-02) 1.309 

P3a 




r(C1-01) 1.337 r(C1-C2) 2.317 
r(C2-01) 1.469 r(C2-02) 1.367 

Figure 4: Calculated structures of pertinent stationary points along 
path C. Distances are given in angstroms (A). 



O- OH \/ 0H 

— I TS4 > I — S 0 TS5 x 0 

^o- — \ 

Int1 Int4 P3a 

Scheme 5: Actual path C obtained by the calculations (as in 
Scheme 3, Int1, TS4, Int4, and TS5 are hydrated by six water mole- 
cules; P3a is the ion-dipole complex with [OH(H20)5]~). 

is shortened, i.e., 1.747 A (TS4), 1.504 A (Int4), and 1.385 A 
(TS5), with a concomitant lengthening of the Cl-01 bond, i.e., 
1.396 A (TS4), 1.429 A (Int4), and 1.531 A (TS5). In the prod- 
uct P3a the proton of the 03-H group is transferred by 
involving the whole water chain to oxygen atom 02. This 
change in the position of the proton is accompanied by 
a shortening and lengthening, respectively, of the Cl-03 
and C2-02 distances: in TS5 r(Cl-03) = 1.383 A and 



r(C2-02) = 1.233 A, while in P3a r(Cl-03) = 1.210 A and 
r(C2-02) = 1.367 A. Finally, it should be noted that intermedi- 
ates of the type Int4 have been proposed [9] to be involved in 
the base-catalyzed reactions of benzocyclobutenediones (reac- 
tion 5 —> 6 in Scheme 2). 

Conclusion 

Ab initio (MP2, SCS-MP2, composite energy approach Eq, 
LPNO-CEPA/1) and density functional methods (M06-2X) 
were applied to study the reaction of cyclobutane-l,2-dione in 
basic solution. The reaction system was modeled by using 
cyclobutane-l,2-dione hydrated with two water molecules and 
[OH(H20)4]~ as the nucleophilic reagent. Three possible reac- 
tion pathways were considered, namely (i) a benzilic acid type 
rearrangement (path A); (b) ring-opening of the bond between 
an aliphatic carbon and that bearing the added OH" group (path 
B); and (c) fission of the bond between the carbonyl carbon and 
that bearing the added OH~ group (path C). Attempts to locate 
path B starting directly from the tetrahedral intermediate Intl 
were unsuccessful. Instead, a reaction sequence diverging from 
the initially formed product Pla of path A was found. Path C 
involved transformations via high energy bicyclic transition 
states and/or intermediates. The final products of these latter 
two paths have comparable (path C) or even substantially more 
negative reaction energies (path B). However, the corres- 
ponding Gibbs free energies of activation are quite large. With 
respect to Intl these are at the CEP A/1 level 22.2 kcal mol -1 
(path B; with respect to Pla AG^ = 44.3 kcal mol" 1 ) and 
44.4 kcal mol -1 (path C). In contrast, path A is not only strongly 
exergonic but also has a significantly lower activation energy, 
AG react = -22.2 kcal mol -1 and AG^ = 16.5 kcal mol -1 with 
respect to Intl. 

Hence, in agreement with experimental observations [5-7], 
cyclobutane-l,2-dione is calculated to react via benzilic acid 
rearrangement to the ring-contracted product, i.e., 1 -hydroxy - 
cyclopropanecarboxylate (2). 

Computational details 

Geometries were optimized by using the M06-2X density func- 
tional [22] and the 6-31+G(d,p) basis set [34,35] and character- 
ized by frequency calculations as minima or transition states. 
For transition states, IRC calculations [36] were also done. 
These geometries were then used for M06-2X, MP2 [23] and 
SCS-MP2 [24] single-point calculations using the 
6-311+G(2df,2p) basis set [37]. For coupled cluster CCSD(T) 
[38] and CEPA-1 [25-27] calculations the 6-31+G(d) and def2- 
QZVPP [39] basis sets, were used, respectively. Initial coordi- 
nates for [OH(H 2 0) 4 ]- were taken from the WATER27 subset 
of the GMTKN30 database [40,41] and reoptimized with M06- 
2X/6-31+G(d,p). [OH(H 2 0) 4 ]- was then placed about 6 A 
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above 1(H20)2 and the combined system again optimized, 
resulting in the ion-dipole complex. An initial structure for TS1 
was obtained from a relaxed (i.e., optimization of all other coor- 
dinates) potential energy scan of the OH-carbonyl-carbon dis- 
tance; this structure was then refined by transition-state optimi- 
zation and further characterized by IRC calculations along both 
directions of the normal mode corresponding to the imaginary 
frequency. The final structures of both IRC calculations were 
then completely optimized. An analogous procedure, usually 
considering several possible reaction coordinates, was used for 
an initial guess of all other transition states. Bulk solvent effects 
(aqueous solution) were obtained by the SMD solvent model 
[21] at the M06-2X/6-3 lG(d) computational procedure. 
Frequencies obtained at the M06-2X/6-3 1+G(d,p) level are 
unsealed. Gibbs free energies are given relative to the 
separated reactants 1(H20)2 and [OH(H20)4]~ and contain a 
1.9 kcal mol -1 correction for the standard state conversion 
1 atm to 1 mol L -1 at T = 298.15 K. Dispersion corrections to 
the M06-2X results were added by Grimme's DFTD3 proce- 
dure [28]. Programs used were ORCA [42], Gaussian 09 [43], 
GAMESS [44], and DFTD3 [45]; MOLDEN [46] 
and MOLEKEL [47] were used for structure building and 
visualization. 

Supporting Information 

Supporting Information File 1 

Detailed computational results and plot of MP2, SCS-MP2 
and M06-2X vs. CEP A AG re i values, pertinent structural 
data, and Cartesian coordinates of all stationary structures. 
[http://www.beilstein-journals.org/bjoc/content/ 
supplementary/1 860-5397-9-64-Sl.pdf] 
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